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We systematically compare filtering methods used to extract topological excitations from lattice 
gauge configurations. We show that there is a strong correlation of the topological charge densities 
obtained by APE and Stout smearing. Furthermore, a first quantitative analysis of quenched and 
dynamical configurations reveals a crucial difference of their topological structure: the topological 
charge density is more fragmented, when dynamical quarks are present. This fact also implies that 
smearing has to be handled with great care, not to destroy these characteristic structures. 



8th Conference Quark Confinement and the Hadron Spectrum 
September 1-6 2008 
Mainz, Germany 

* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos. sissa. it/ 



Comparison of filtering methods 



F. Gruber 



1. Filtering methods 

Many methods have been developed to extract the IR content from lattice data. Unfortunately, 
all these methods introduce ambiguities and parameters. Thus, to get a coherent picture of the topo- 
logical structure of the QCD vacuum, it is necessary to find ways of controlling or even removing 
these ambiguities. 

One of the first attempts to filter out the UV "noise" has been APE smearing [1], defined as: 
^ APE = P su(N c ) {(1 " a APE )Uf + ^(staples)} , (1.1) 

where (Xape determines the weight of the old link and the sum of the attached staples. The right 
hand side has to be projected back to the gauge group. Unfortunately, there is no unique mapping. 
One approach is to take the unitary part of the polar decomposition and normalize this matrix by 
its determinant. Stout smearing [2] circumvents this projection by using the exponential map: 

= exp{^(t/,p MV )} • Uf, (1.2) 

where Qii{U ,p^ v ) is a hermitian matrix constructed from all plaquettes containing the old link 
and weighted by factors p^ v . We use the common choice p^ v = pstout for isotropic smearing. 

A relatively new method is Laplace filtering [3]. The filtered links are obtained from a spectral 
sum of the lowest eigenmodes of the covariant lattice Laplacian 1 : 

U^ lace (x) =P su{Nc) j- £ • (1-3) 

This procedure acts as a low-pass filter in the sense of a Fourier decomposition. At this point it 
should be stressed that Laplace filtering is completely different from smearing, because it is based 
on rather global objects, namely the eigenmodes, and does not locally modify the gauge links in 
contrast to smearing. 

Taking the filtered links as a starting point, one can reconstruct the topological charge density 
q(x) = Tr {F tlv {x)F ilv {x)) / I6n 2 from an improved field strength tensor [4]. 

Also the fermionic definition of the topological charge, via the eigenmodes y of a chiral 
Dirac operator, has been used to explore the IR structure [5]. For this so called Dirac filtering one 
truncates the sum in 

N / 3 \ 

q Dirac (x) = £ f — - 1 J 1/4 Wn (*) (1.4) 

and takes only the lowest N modes into account. While the zero-modes determine the total topo- 
logical charge Q = Y,q{ x ) due to the index theorem, the non zero-modes modify the local structure 
of the density, leaving the total charge unaffected. 

2. Comparison of the different methods 

In an earlier study a qualitative and quantitative similarity of the introduced filtering meth- 
ods for quenched SU(2) gauge configurations has been observed [6]. One central element of this 

'The original link is reproduced for all eigenmodes, N = N c - Vol, with no projection needed. 
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comparison is the correlator of two topological charge densities qA(x) and q B (x) defined by: 

Xab = (1/V)£ (qB(x)-q B ), (2.1) 

where the mean values are subtracted for convenience. From this we can construct a quantity that 
reflects the "matching" of two methods: 

Sab = (2.2) 
Xaa Xbb 

Sab is obviously equal to one, if qA(x) is proportional to q B (x) and deviates the more from one, the 
more the densities differ. 

The main idea is now to relate different filter parameters for those combinations where S is 
maximal. In fig. 1 the contour lines of S for several methods and parameter ranges are shown. On 
the right hand plot two exemplary combinations are indicated that correspond to the best matching 
value for different filtering strengths. 

An interesting observation is that there is an almost one-to-one correspondence for n steps of 
APE and n steps of Stout smearing when GCape ~ 6 • Pstout- As seen in the plot on the lhs. of fig. 1, 
S > 0.95 for a large number of smearing steps. This is consistent with results by Capitani et al. [7], 
where such a relation has been derived from perturbation theory. While they have focused on global 
observables with up to 3 smearing steps, our nonperturbative result reflects the local similarity of 
both methods and their strongly correlated topological charge densities up to 50 steps. 
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Figure 1: Level curves of S = 0.95,0.85, . . . (starting from the diagonal) for APE vs. Stout smearing (left) 
and E = 0.8,0.7, . . . (starting from the inside) for APE smearing (S) vs. Laplace modes (L) (right). ▲ and • 
mark two examples of "matching" parameters for weak and strong filtering respectively (from [6]). 



3. Cluster analysis of the topological charge density 

Another important challenge is to extract observables from lattice data, that could be com- 
pared with continuum models of the vacuum. One possibility is to analyze the cluster structure of 
the topological charge density. Two lattice points belong to the same cluster, if they are nearest 



3 



Comparison of filtering methods 



F. Gruber 



0.9 • 



0.8 



? 0.7 



'0.6 



0.5 



0.4 



dynamical 
quenched 



instanton gas, N f -2 



instanton gas, quen. 



dynamical • 
quenched 



io J | 

CD 



O. 

c 



10" 



8 9 10 15 22 

smearing steps 



30 



40 



6 7 



8 9 10 15 22 

smearing steps 



30 



40 



Figure 2: left: Exponent t, of the analysis for clusters common to APE and Stout smearing. The solid lines 
show the values predicted from the dilute instanton gas. right: Total number of distinct clusters for a constant 
fraction / = 0.0755 of points lying above the cut-off. Less than 6 steps are not considered, as the definition 
of the topological charge gets ill-defined. Errors have been calculated using an ensemble average over 10 
configurations but are partly too small to see. 



neighbors and have the same sign of the topological charge density. Bruckmann et al. [6] found 
a power law for the number of clusters as function of the ratio of points with \q{x)\ lying above 
a variable cut-off q cut and the total number of lattice points. The exponent <?; of this power law is 
highly characteristic for the topological structure of the QCD vacuum. Different models lead to 
different predictions, which allows for a very sensitive test. If one has for instance pure noise, the 
exponent is 1, as every point forms its own cluster. On the other hand one will have an exponent 
close to zero for very smooth densities with large structures. 

To reduce ambiguities we take only those clusters into account, which are common to different 
niters, whose parameters were matched according to maximal values of S. So, if there is an artifact 
coming from one method, it is unlikely that this artifact will also be seen by the other. 

The exponent for clusters common to APE and Stout can be found in fig. 2 (left). We used one 
quenched and one dynamical Nf = 2 ensemble with equal lattice spacing (see tab. 1). Obviously 
the exponents of the dynamical configurations lie above the quenched values. 

In order to interpret the cluster exponent, a model of dilute quantized topological objects of 
general shape and with a size distribution d(p) ~ has been considered in [6]. It leads to E, = 
(1 +4/Q3 + (in 4 dimensions). Following this model, our findings give a larger coefficient j8 
in the dynamical case. Hence, smaller topological objects become suppressed. 

Moreover, the rhs. of fig. 2 shows that for a fixed number of points, lying above the cut-off, 
much more clusters are found in the dynamical case. Thus we conclude that when fermion loops 
are taken into account, the topological structure is more complex and fragmented, in the sense of 
larger number of distinct objects per volume. This seems to be in accordance to the findings of 
the Adelaide group, where small instantons have been seen to be suppressed in the presence of 
dynamical quarks, while the total number of instantons increased, see fig. 6 in [8]. 

The difference of the cluster exponents quenched vs. dynamical vanishes for stronger smearing 
(~ 30 steps) and the exponents settle down to the same plateau. So we have reasons to believe that 
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too much smearing destroys the impact of dynamical quarks. 

On the lhs. of fig. 2 we have included for comparison the exponents t, =7/11 « 0.64 and 
% = 23/35 w 0.66 for the SU(3) instanton gas without resp. with dynamical quarks. Taking the 
dilute instanton gas as a simplified model, it is obvious that the true vacuum should have a higher 
exponent, as more structures are present. However, the result in fig. 2 (left) shows that this is only 
the case for very few smearing steps, for slightly stronger filtering we reach smoother configura- 
tions than predicted by the dilute instanton gas. This is another indication of smearing artefacts. 





lat. size 


lat. spacing 


Plw 


m 


quenched 


16 3 -32 


0.148 


7.90 




dynamical 


16 3 -32 


0.150 


4.65 


-0.060 



Table 1: Ensembles were generated with the Luscher-Weisz gauge action and a chirally improved Dirac 
operator [9]. For the dynamical simulations two flavors of mass degenerate light quarks were used [10]. 

4. Conclusion and outlook 

In conclusion, we have found a strong correlation of the topological charge densities obtained 
from APE and Stout smearing. Furthermore, our first results for dynamical quarks imply that the 
topological structure is more complex and fragmented in the presence of fermion loops. But there 
are also indications that smearing has to be used with great caution, especially when dealing with 
dynamical configurations. The smallness of the cluster exponent is a sign that smearing is more 
destructive in SU(3) than in SU(2). Preliminary results for clusters common to APE smearing and 
Laplace filtering do not show such artefacts, as Laplace filtering preserves smaller objects better. 
This effect is under investigation. 
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